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ABSTRACT 

In blazars snch as 3C 279, GeV gamma-rays are thonght to be prodnced by 
inverse Compton scattering of soft photons injected from external sonrces into 
the jet. Becanse of the large bnlk Lorentz factor of the jet, the energy of soft 
photons is Doppler shifted in the comoving frame of the jet, and the scattering is 
likely to occnr in the Klein-Nishina regime. Althongh the Klein-Nishina effects 
are well known, the properties of the electron and emission spectra have not been 
stndied in detail in the environment of blazars. We solve the kinetic eqnation 
of electrons with the spatial escape term of the electrons to obtain the electron 
energy spectrnm in the jet and calcnlated the observed emission spectrnm. In 
calcnlations of the Compton losses in the Klein-Nishina regime, we nse the dis¬ 
crete loss formalism to take into acconnt the signihcant energy loss in a single 
scattering. Althongh the scattering cross section decreases becanse of the Klein- 
Nishina effects, ample gamma rays are emitted by inverse Compton scattering. 
When the injection spectrnm of electrons obeys a power law, the electron spec¬ 
trnm does not follow a broken power law, as a resnlt of the Klein-Nishina effects, 
and a large nnmber of high-energy electrons remain in the emitting region. 

Subject headings: gamma rays: theory — radiation mechanisms: nonthermal 


1. Introduction 


Inverse Compton (IC) scattering is an important process in prodncing high-energy pho¬ 
tons in astrophysical objects such as active galactic nuclei and gamma-ray bursts. Blazars 
are one of the important sites at which IC scattering is a major radiation process. The radia¬ 
tion from blazars is well described by synchrotron emission and IC scattering (e.g., Sikora et 
ah 1994; Inoue & Takahara 1996; Ghisellini et al. 1996). The spectral energy distribution 
of blazars has two peaks in the representation; one is in the range between IR and 

X-rays, and the other, in gamma rays such as GeV and TeV regions. The emission region of 
blazars is thought to be in the relativistic jet close to the central black hole, and the nonther- 
mal electrons accelerated in the jet are probably responsible for the emission. Synchrotron 
radiation emitted by the nonthermal electrons produces the low-energy peak, and IC scat¬ 
tering by the same nonthermal electrons produces the high-energy peak. Although there 
are hadronic models that assume relativistic protons (e.g., Mucke et al. 2003), we conhne 
ourselves to leptonic models, which assume radiation from nonthermal electrons/positrons 
alone. 

In the jet of 3C 279, nonthermal electrons Compton scatter external soft photons to 
produce GeV gamma rays (Inoue & Takahara 1996); the source of the external soft photons 
may be accretion disks or broad-line regions (e.g., Dermer & Schlickeiser 1993; Sikora et al. 
1994). In a previous paper (Kusunose et al. 2003), we showed that the electron spectrum 
in the emission region of 3C 279 does not obey a conventional broken power law, which 
appears if scattering occurs in the Thomson regime. The electron energy spectrum of 3C 
279 is found to be harder than the broken power law in the high-energy region because of 
the Klein-Nishina (KN) effects. 

When the electron Lorentz factor 7 and the seed soft photon energy normalized by 
the electron rest mass x = hv/{meC^) satisfy 7 a; < 1, IC scattering occurs in the Thomson 
regime. Here me and c are the electron rest mass and the speed of light, respectively. On 
the other hand, if 7 a; > 1, IC scattering occurs in the KN regime (e.g., Blumenthal & 
Gould 1970; Rybicki & Lightman 1979). In astrophysical objects high-energy radiation 
is often emitted by nonthermal electrons with a power-law distribution. In the objects in 
which nonthermal electrons are injected, radiate, and escape, the electron spectrum becomes 
steeper for 7 > 7 br, where 7 br is the Lorentz factor determined by the balance between the 
Compton cooling and the electron escape, if the scattering occurs in the Thomson regime. 
When the injection spectrum obeys a power law, qe{'^) oc 7 “^, the steady state spectrum 
obeys a broken power law, ne( 7 ) oc 7 “^ for 7 < 7 br and ne( 7 ) oc 7 “^*“^ for 7 > 7 br- 

The effects of IC scattering in the KN regime have been studied by several authors. 
Blumenthal (1971), for example, showed that the electron spectrum becomes harder in 
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the KN regime, assuming that monoenergetic soft photons are scattered. Zdziarski & Krolik 
(1993) demonstrated that when the electrons are injected with a power-law distribution with 
the power-law index of 2, the emission spectrum is flat in the representation if the KN 

effects are taken into account, although they assumed that the scattering cross section is zero 
for '-yx > 3/4 to simulate the KN effects. These properties were applied to the explanation of 
the TeV emission from Mrk 421 and its flat spectrum in the GeV—TeV region in the v-vFy 
representation. Zdziarski (1989) solved kinetic equations for various types of the injection 
spectra of electrons and soft photons, using the KN cross section. 

In above work (Blumenthal 1971; Zdziarski 1988), the electron escape from the emis¬ 
sion region was not included. In objects such as blazars, the emission region exists in the 
relativistically moving jets. Because the advection or adiabatic expansion of nonthermal par¬ 
ticles removes the particles from the emission region on a timescale close to the light crossing 
time i?/c, where R is the radius of the emission region, the effect of escaping particles on 
the spectra of electrons and photons is important. In this paper we show the properties of 
the spectra of electrons and photons when the electron escape is included and the scattering 
occurs dominantly in the KN regime. We also assume that synchrotron radiation is not im¬ 
portant for radiative cooling; this is justified when ample soft photons from external sources 
are injected as seed photons of IC scattering. 

In §2 the kinetic equation of electrons is presented, and numerical results are given in 
§3. Finally we summarize our results in §4. 


2. Kinetic Equation 

The kinetic equation for the steady state electrons is given by 

-ne(7) [ di C(j, Y) + f dYUeiY) C{Y, 7 ) + Qeil) - = 0 , ( 1 ) 

Jl J 'y ^esc 

where nYY is the number density of electrons per unit 7 , qeY) is the injection rate of 
nonthermal electrons per unit volume and 7 , tesc is the electron escape time, and C(j,y) 
describes the transition rate of an electron from 7 to 7 '. We assume that geij) is given by 

QeY) = go 7“^ exp(-7/7„iax) , 7 > 7min , (2) 

where p, 7 min, 7 max, and go are parameters. 

In scattering the electron Lorentz factor 7 decreases to 7 '. When the number spectrum 
of soft photons per unit volume, nsoft(x), is given, the rate of scattering (7 — 7 ') per electron 
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is given by Jones (1968; see also Zdziarski 1988) as 

nsoft{x) dx, 

where 

E^^x, = E>E., r = i(l + 4j, (4) 

Here ut is the Thomson cross section, and we assume that the soft photons are isotropic. 
Equation (3) retains the zeroth-order term in a double expansion of the exact rate in l/y^ 
and x/'j] the exact cross section is given by Jones (1968; see also Coppi & Blandford 1990). 
Since we are interested in the KN regime, the above approximated cross section is enough 
to evaluate the electron spectrum for 7 3> 1. We assume that soft photons with energy Xs 
are isotropically distributed around the jet, where Xs is measured in the frame of the soft 
photon source. In the jet frame, the soft photon energy is Doppler shifted to a; = Tx^, where 
r is the bulk Lorentz factor of the jet. The radiation held is approximated to be isotropic 
in the jet frame for simplicity. 


(7(7,7) = 


3(Jt c 
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To solve equation (1), we discretize the equation as 


i—1 J 

C'(7i,7i)7iAln7+ ^ ne(7i)C'(7j,7j)7jAInq + ^^(7^) 
j=l j=i+l 


ne{li) 

fesc 


0 , 


( 5 ) 


where 7 * {i = 1, • ■ ■ , J) is binned logarithmically and Alny = In 7^+1 — ln 7 j with A In 7 
being constant (typically Alny = 2 x 10“^). Here we omitted the term ( 7 ( 7 *, 7 *), because 
most electrons lose energy in scattering. We then obtain 


ne(7i) = A ^ 


■ J 

^ ne(7j)C(jj, ji)jj Ain J + qe(ji) 

J=i+1 


( 6 ) 


where 

1 

A = Y^ (7(7*, 7i)7i A In 7 + — . (7) 

j=i 

We calculate Uei'-fi) numerically from equation ( 6 ) with the upper boundary condition such 
that Uei'-fj) = 0. Because ( 7 ( 7 , 7 ') equation (3) is applicable to large values of 7 , our 
numerical calculations are limited to 71 > several. 


The injected electrons lose energy by IC scattering and move to the lower energy region 
in the energy space. The cooling rate of IC scattering is given by 

7 = ^ c?7'(7'-7)C'(7,7)- 


( 8 ) 
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The cooling time is then dehned as tic( 7 ) = 7 /l 7 l- the Thomson regime, tic( 7 ) decreases 
with increasing 7 and attains the minimum value at 7 ~ Xq where Xq is the characteristic 
energy of soft photons. For 7 > Xg where the scattering occurs in the KN regime, tic ( 7 ) 
increases. 


Nonthermal particles in the jet may escape from the emission region by advection or 
diffusion, or adiabatic expansion effectively removes them. Because of the lack of knowledge 
on the electron escape time, tesc is assumed to be independent of 7 and given by tesc = 
Cesc iiiiii{tic( 7 )}, where Cesc is a parameter and min{tic( 7 )} is the minimum value of tic ( 7 ), 
which is taken at 7 ~ Xq^. 


For comparison with the scattering in the KN regime, we also calculate the electron 
spectrum derived by the Thomson approximation. The kinetic equation of electrons in the 
Thomson regime is given by 

;^[7t( 7) «e(7)] = 9e(7) - , (9) 

«7 tesc 

where the cooling rate is calculated (Lightman & Zdziarski 1987) 

/4 \ r 3/(47) 

7 t ( 7 ) = - o-tc ( - 7 ^ - 1 ) j nsoit{x)xdx . (10) 

Here the scattered photon energy x' is assumed to be given by x' = ( 4 / 3 ) 7 ^a;, and the 
scattering cross section is assumed to be zero for xy > 3/4 to simulate the scattering in the 
KN regime. The solution of equation (9) is given by 


«e(7) 


7 t ( 7 ) 


Qeii) exp 


dr] 

ixiv) 




( 11 ) 


3. Numerical Results 

We solve equation (1) numerically to obtain the electron spectrum and calculate the 
emission spectrum. The emission region is assumed to be in a relativistically moving blob, 
and we calculate the emission spectrum according to Georganopoulos et al. (2001), who 
includes the effects of the beaming and the KN cross section, because our interest is in the 
application of our model to blazars. We hrst solve equation (1) without the term of electron 
escape to compare with the results with the escape term. In Figure 1 we show the electron 
spectra for different values of 7 ma v and p = 2 and 1.8 (p = 1.8 was used in our model of 3C 
279; Kusunose et al. 2003). In this hgure, nsoft(a^) = noS{x — Xq) and XorrieC^ = 1.25keV, 
where S is the Dirac delta function; when the bulk Lorentz factor of the jet is F = 25, 
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this soft photon energy is 50 eV in the frame of the soft photon source. The value of go 
is chosen to set the injection rate of electrons to 1 cm“^s“^, and the value of no is hxed 
to set the energy density of soft photons to 8.3 x 10“^ ergs cm“^; this value of the energy 
density corresponds to 1.3 x 10“^ erg cm“^ in the frame of soft photon source if T = 25. 
When p = 2, ne{'y) oc 7 “^ for 7 x 0 < 0.1, i.e., scattering occurs in the Thomson regime. The 
KN effects appear for 7 x 0 > 0.1, and the electron spectra become harder. For example, 
when 7 max = 10 ®, ne( 7 ) oc 7 “^-^ for p = 2 and 10 < 7 x 0 < 10 ^, with a prominent peak at 
7x0 ~ 3 X 10^. It should be noted that if we solve equation (1) at the lower energy region 
down to 7 = 1 , the electron spectrum should have an inhnite peak at 7 = 1 because of 
cooling. Since we solved equation ( 1 ) at large values of 7, this peak does not appear in the 
hgures. 

The deviation of the Thomson approximation (eq. [11]) from the KN formulation (eq. 
[1]) is shown in Figure 2. Here the injected soft photons obey the Planck distribution with a 
temperature of 1.25 keV and an energy density of 8.3 x 10“^ ergs cm“^. The values of ne(7) 
are different from those in Figure 1, because the external soft photon spectra are different. 
Since electrons with 7 > 3/(4a;) do not suffer from cooling in the Thomson approximation, a 
large number of electrons stay in the high-energy region. Thus, just assuming that the cross 
section vanishes for 7 > 3/(4a;) considerably overestimates the number density of electrons 
in the high-energy region. More accurate treatments are needed to properly estimate the 
high-energy electron spectrum. 

The observed photon spectra emitted by the electrons with p = 2 shown in Figure 1 are 
presented in Figure 3 (the spectra from electrons with p = 1.8 given in Fig. 1 are shown in 
Fig. 4). It is assumed that the emission region is in the moving blob with F = 25 at redshift 
= 0.538 (the redshift of 3C 279) and that the size of the emission region is i? = 7 x 10^^ cm 
(the blob size appropriate for 3C 279 in our model; Kusunose et al. 2003). The soft photon 
spectrum of the external photon source is monochromatic, with energy XQUieC^ = 1.25 keV. 
Although the scattering cross section decreases owing to the KN effects for 7 > O.lxg^ 
there are still ample hard photons to be observed. In these hgures synchrotron spectra with 
B = 0.3 G are also shown. The energy density of the magnetic held is 3.6 x 10“^ ergs cm“^, 
which is to be compared with the soft photon energy density of 8.3 x 10“^ ergs cm“^. It is 
clear that when 7max increases, the synchrotron luminosity also increases. On the other hand, 
the Compton luminosity does not, because of inefficient cooling in the KN regime. Thus if 
the Compton cooling occurs mainly in the KN regime, the dominance of the synchrotron 
luminosity over the Compton luminosity does not mean that the magnetic energy density 
dominates other energy contents such as external soft photons. It should be noted that 
the cooling rate of synchrotron emission is not taken into account in the electron kinetic 
equation and that the synchrotron spectra are shown only for the purpose of seeing the 
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effect of different shapes of the electron spectra in the high-energy region. For 7 max ^ 10^, 
the synchrotron cooling dominates the Compton cooling when the magnetic held is hxed at 
0.3 G. On the other hand, if i? ~ 0.02 G, the synchrotron cooling time is one tenth of the 
Compton cooling time for 7 < 10®. 

As mentioned above, we solved equation ( 1 ) at large values of 7 . If we solve the equation 
at small values of 7 down to 1 , there should be inhnite number of electrons at 7 = 1 , and 
they should produce a peak at hv r^hz/Q 31keV in the photon spectrum, where hvQ is 
the soft photon energy in the rest frame of the soft photon source: 50 eV in our numerical 
calculations. 

From Figures 3 and 4 it is found that the small differences in the electron spectrum 
caused by the different values of 7 max result in different values of the power-law index of 
the lower energy part of the synchrotron spectrum, or the part of the spectrum below the 
synchrotron peak energy. This might be used to determine observationally the KN effects 
on the electron spectrum by the observations in radio through optical bands. 

The electron spectrum in the jet with electron escape from the emission region is shown 
in Figure 5 for the different values of Cesc, assuming that the soft photons are monochromatic. 
In this hgure, 7 max = 10"^, p = 2, and nieC^xo = 1.25 keV are assumed. We assumed that 
Cesc oc tesc, and Cesc = 1 correspouds to tesc = 3.6i?/c for the adopted parameter values. When 
the escape time is shorter, the electron density decreases, the spectrum becomes flatter, and 
the peak energy of the bump decreases. When Cesc = 0.1, the shape of the electron spectrum 
is almost the same as that of the injected electrons. The electron spectra for the different 
values of 7 max are shown in Figure 6 with Cesc = 1- Even when the electron escape time is 
a few times the light crossing time {R/c), the energy density of high-energy electrons with 
7 > Xq ^ is still large, and the emission from those electrons are expected to be as shown 
below. 

The electron spectra with the Thomson approximation and those with the KN cross 
section are compared in Figure 7. The injected soft photons follow a blackbody distribution 
with a temperature of kT = 1.25 keV and an energy density of 8.3 x 10“^ ergs cm“®. As with 
the case in which there is no electron escape, the Thomson approximation results in more 
electrons in the high-energy region than the calculations with the KN cross section. In the 
Thomson approximation there is a dip at 7 ~ 100, which corresponds to 'ykT / {rUeC^) ~ 0.24. 

The photon spectra from the jet with electron escape are shown in Figures 8 and 9. 
In both hgures, the external soft photons are monochromatic. Synchrotron emission with 
B = 0.3 G is also shown for comparison with the IG emission. In Figure 8 , 7 max = 10^, 
p = 2, and various values of (esc sue assumed. When (esc = 0.1, the Gompton luminosity 



is higher than the synchrotron Inminosity. When (esc increases, both the Compton and 
synchrotron Inminosities increase, while the spectral shape of both components changes as 
a resnlt of the change in the electron energy distribntion dne to the KN effects. In Figure 
9, Cesc = 1 ; P = 2 , and the different values of are assumed. It is found that the peak 
of uFiy of the Compton spectrum occurs at hu ~ lO^meC^ for large enough values of 7 max 
(> 10"^) when Cesc = 1- When Cesc is larger, the peak of uFi, occurs at a slightly larger value 
of photon energy. The peak energy hu ~ lO^mgC^ is due to the electrons with 7 ~ 160, 
i.e., 7 x 0 ~ 0.4, which scatter seed photons in the Thomson regime producing photons with 
energy ~ Taking into account = 0.538 and T = 25, the observed peak energy 

appears at hu ~ lO^mgC^. Because the synchrotron luminosity is larger than the Compton 
luminosity for 7 ma x > a few x 10 *^, the properties of electron and photon spectra described 
here for 7 max ^ a few x 10^ are correct when B is smaller than 0.3 G to keep the synchrotron 
luminosity below the Compton luminosity. 


4. Summary 

We studied the effects of inverse Compton (IC) scattering in the Klein-Nishina (KN) 
regime, paying attention to the application to the gamma-ray emission from blazars. For 
this purpose we solved the kinetic equation of electrons with electron escape, assuming that 
electron cooling is predominantly by IC scattering of external soft photons. The photon 
spectrum was calculated, assuming that the emitting region is moving relativistically along 
the line of sight, i.e., emission from a relativistic jet is considered. The seed soft photons of 
IC scattering are assumed to come from external sources. Because of the Doppler effect, the 
energy of the seed photons is higher for electrons in the jet. Then the KN effects appear for 
smaller values of 7 compared with electrons in a stationary emitting region. In our numerical 
calculations, the bulk Lorentz factor is set to be 25, and the seed soft photon energy is 50 
eV in the frame of the soft photon source, which is based on our model of 3C 279 (Kusunose 
et ah 2003). For these parameters, electrons with 7 > 400 scatter soft photons in the KN 
regime. 

When the escape of electrons can be neglected compared with the cooling, the electron 
spectrum becomes harder than the injection spectrum for 7 > Xq^, as Blumenthal (1971) 
showed, because of the less efficient cooling in the KN regime. Although the cooling rate is 
smaller compared with that in the Thomson regime, gamma rays with energy of ~ 7 max are 
amply emitted by IC scattering. When electron escape is taken into account, the number 
of high-energy electrons decreases. However, gamma rays in GeV region are still emitted by 
IC scattering. 



9 


In astrophysical objects such as blazars, where IC scattering is a major radiative pro¬ 
cess, the spectrum of nonthermal electrons needs to be carefully considered. Although a 
simple power-law spectrum of electrons might £t the observed emission spectrum, it may 
not represent the physical conditions in the source. The self-consistent treatment of photons 
and particles should be taken into account. We hnally remark on the effects of synchrotron 
radiation and synchrotron-self-Compton scattering on the electron spectrum. When the 
magnetic held is strong enough, these processes decrease the number density of electrons 
with very large Lorentz factors, and the upturn in the high-energy tail of the electron spec¬ 
trum becomes less prominent. This effect was included in our previous paper (Kusunose et 
ah 2003), where the KN effects are still apparent in the electron spectrum. 


This work has been partially supported by Scientific Research Grants (M.K.:15037210; 
F.T.: 14079205 and 16540215) from the Ministry of Education, Culture, Sports, Science and 
Technology of Japan. 
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Fig. 1.— Electron spectra for various values of 7max without particle escape from the emission 
region. Thick and thin lines are for p = 2 and 1.8, respectively. The injection of monochro¬ 
matic soft photons with rrieC^XQ = 1.25keV is assumed. The KN effects are apparent for 
72^0 > 0.1. 
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Fig. 2.— Electron spectra calculated with the KN cross section {thick lines) and with the 
Thomson approximation {thin lines). Injected electrons obey a power law with p = 2, and 
the external soft photons follow a blackbody distribution with a temperature of 1.25 keV. 
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Fig. 3.— Photon spectra emitted by electrons with p = 2 shown in Fig. 1 and various 
values of 7max- The spectra are shown in the observer’s frame. The curves on the right are 
the emission by IC scattering and the curves on the left are by synchrotron radiation with 
B = 0.3 G. 
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Fig. 5.— Electron spectra with various values of escape time; 7max = 10"^ and p = 2 are 
assumed for the injection spectrum of electrons. The soft photons are monochromatic with 
the energy 1.25 keV. From upper to lower curves, Cesc = oo (no escape), 10, 5, 3, 2, 1, and 

0 . 1 . 
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Fig. 7.— Electron spectra with escape calculated with the KN cross section {thick lines) 
are compared with those calculated by the Thomson approximation {thin lines). Electron 
escape time is tesc = 1-5 x 10^ s, which corresponds to Cesc = 1 for the KN case. 
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Fig. 8.— Photon spectra emitted by the electrons shown in Fig. 5. Curves on the left are 
by synchrotron radiation with B = 0.3 G, and curves on the right are by IC scattering. 









Fig. 9.— Photon spectra emitted by the electrons shown in Fig. 6 (Cesc = 1) produced by 
IC scattering and synchrotron radiation with B = 0.3 G. 











